Driven Disordered Periodic Media with an Underlying Structural Phase Transition 
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We investigate the driven states of a two-dimensional crystal whose ground state can be tuned 
through a square-triangular transition. The depinning of such a system from a quenched random 
background potential occurs via a complex sequence of dynamical states, which include plastic 
flow states, hexatics, dynamically stabilized triangle and square phases and intermediate regimes of 
phase coexistence. These results are relevant to transport experiments in the mixed phase of several 
superconductors which exhibit such structural transitions as well as to driven colloidal systems 
whose interactions can be tuned via surface modifications. 
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The description of phases and phase transitions in 
driven steady states is a central theme in the statis- 
tical mechanics of non-equilibrium systems 0. A vari- 
ety of such states are obtained in the depinning and 
flow of randomly pinned, periodic media, such as charge- 
density wave systems and Abrikosov flux-line lattices in 
the mixed state of type-II superconductors |2j. Many 
such superconductors exhibit structural transitions in the 
mixed state, typically between flux-line lattices with tri- 
angular and rectangular symmetry [j| . Colloidal systems 
such as PMMA spheres coated with a low-molecular- 
weight polymer undergo a remarkable variety of solid- 
solid transformations in an external field^]. While apply- 
ing a sufficiently large current depins the flux lines from 
the quenched random disorder present in all real mate- 
rials, the possibility of driving colloidal particles in two 
dimensions across a disordered substrate has also been 
raised 5] . What links these diverse systems is the generic 
problem of understanding the competition between an 
underlying structural phase transition in a pure periodic 
system as modified by disorder, and the non-equilibrium 
effects of an external drive. This Letter proposes and 
studies a simple model which describes this physics. 

Our model system is two-dimensional and consists of 
particles with two and three-body interactions 6J . The 
three-body interaction, parametrized through a single pa- 
rameter U3, tunes the system across a square-triangular 
phase transition. Our central result, the sequence of 
steady states obtained as a function of increasing force 
F for various values of V3, is summarized in the dynam- 
ical "phase" diagram of Fig. We obtain a variety of 
phases: pinned states which may have dominantly trian- 
gular or square correlations, a moving liquid/glass phase, 
a moving hexatic glassy phase, flowing triangular and 
square states ordered over the size of our simulation cell 
and a dynamic coexistence regime between these ordered 
phases. We discuss our characterization of these states 



and the applicability of simple dynamical criteria for non- 
equilibrium phase transitions between them. 

The model: Particles interact in two dimensions 
through the interaction potential 1/2 Y^ijtj ^M^y) + 



Vs(ri, rj, 7*fc), where is the position vector 
= \ rij \ = \Tj-Til V 2 (rij) = v 2 (^) 12 and 
V3{fijSin 2 (46 ijk )fj k + permutations] 0. 
/) = [rij - r ) 2 for r. tj < 1.8a 



1/6 E 
of particle i, 

V 3 (ri,rj,rk) = 
The function /y = /(fy, 
and otherwise and 9ijk is the angle between Tji and 
Yjk- The two-body (three-body) interaction favours a 
triangular (square) ground state. Energy and length 
scales are set using v 2 = 1 and er = 1. Particles 
also interact with a quenched random background mod- 
eled as a Gaussian random potential^ V<j(r) with zero 
mean and exponentially decaying (short-range) corre- 
lations, The disorder variance is set to in = 1 and 
its spatial correlation length is £ = 0.12. The system 
evolves through standard Langevin dynamics; = Vj 
and Vj = fi' nt — avi + F + T]i(t). Here v, is the veloc- 
ity, fi' nt the total interaction force, and r]i(t) the random 
force acting on particle i. A constant force F = {F X ,Q} 
drives the system and the zero-mean thermal noise r\i (t) is 
specified by < ^(^(t') >= 2TS ij 8(t-t') with T = 0.1, 
well below the equilibrium melting temperature of the 
system. The unit of time r = a<7o/i>2, with a — 1 the 
viscosity. 

Simulation details: Our system consists of 1600 par- 
ticles in a square box at number density p = 1.1. At 
T = 0.1, the pure system remains triangular upto W3 = 
1.5. For larger a square phase is obtained. Larkin 
length estimates|2j yield L a /a ~ 100, with a = l/p 1 ^ 2 
the lattice parameter, somewhat larger than our system 
size. 

Equilibration in the disorder potential is achieved 
through a simulated annealing procedure in which the 
disorder potential is varied by increasing its strength 
from zero to unity in small steps. The system is stabi- 
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FIG. 1: The dynamical phase diagram in the v$ — F plane. 
The phases are: A - pinned triangle, B - pinned square, C 
- plastic flow/isotropic liquid, D - moving hexatic, E - mov- 
ing triangle, F - moving square and G - dynamical square- 
triangle coexistence. Two points with error bars show the 
boundary of G for = 6, obtained by averaging over 24 
disorder realizations, as an example. The boundaries for all 
other transitions are considerably sharper. The inset shows 
the centre of mass velocity « cm as a function of the driving 
force F as the force is increased (bold) and then decreased 
across the depinning transition. The right panel shows S(q) 
for the plastic flow state (i), hexatic glass (ii), triangular solid 
(iii) and square solid phases (iv) at V3 = 6.0 and F = 10, 20, 60 
and 140 respectively. To obtain S(q), 50 independent config- 
urations were used. The structure in (iv) reflects the presence 
of two mutually misoriented square crystallites. 

lized for a minimum of 2 x 10 5 Monte Carlo steps at each 
disorder strength. Such annealed configurations are our 
initial inputs to the Langevin simulations. We evolve the 
system using a time step of 10 _4 r. The external force P 
is ramped up from a starting value of 0, with the system 
maintained at upto 10 s steps at each P. 

Observables: We monitor structural observables, such 
as the static structure factor S(q) = J2ij exp(— iq • r^). 
Delaunay triangulations yield the probability distribu- 
tions P(n) of n = 4,5,6 and 7 coordinated parti- 
cles i^}2 n P{n) = 1)8]. We define order parameters 
ijj = (P(4) - P(6))/(P(4) + P(6)) to distinguish be- 
tween square and triangular phases, ipt\ — (P(6) — 
P(5)-P(7))/(P(6) + P(5) + P(7)) to distinguish between 
liquid (disordered) and triangular crystals and tpn — 
(P(4) - P(5) - P(7))/(P(4) + P(5) + P(7)) to distinguish 
between liquid and square crystals. In addition, we com- 
pute the hexatic order parameter ips = Y^ij exp(— iQOij) 
and its correlations, where 9 is the bond angle measured 
with respect to an arbitrary external axis. The dynami- 
cal variables we study include the center of mass velocity 
v cm , the particle flux and its statistics and the Koshelev- 
Vinokur (KV) "shaking temperature" T" Q] appropriate 
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FIG. 2: (a) The shaking temperatures for a system with 
vz — 6 as a function of F in the drive (x) and transverse 
(y) directions. These are obtained by averaging over 100 in- 
dependent configurations as well as over 25 separate disorder 
realizations. For F in the coexistence region there is a signifi- 
cant enhancement of the shaking temperature in excess of the 
prediction of KV (Ref. 9]). The lines represent fits to the KV 
form, (b) The hexatic correlation function ga(r) at varying 
force values at 173 = 6, averaged over 50 independent configu- 
rations, (c) P(n) for n = 4, 5, 6, 7 (see text) vs F averaged as 
in (a). 




FIG. 3: (color-online) A single particle configuration show- 
ing square-triangle coexistence at F = 107 together with the 
computed Delaunay mesh. The particles are colored accord- 
ing to the number of neighbors n = 4 (magenta), 5 (green), 
6 (blue) and 7 (orange). Particles with coordination 5 are 
present mainly in the interfacial region while those with 7 are 
associated with isolated dislocations. 



to the drive and transverse directions and obtained from 
T? = {[v<*-v2 m f)/2,a = x,y. 

For small F the solid is pinned. A disorder broadened 
version of the equilibrium triangle (A) to square (B) tran- 
sition results as V3 is varied across the T = transition 
value at small F; here and below alphabets in brackets 
refer to the states labelled in Fig. ^ The triangular (A) 
phase is favoured at non-zero F. Upon further increas- 
ing P, the system undergoes a discontinuous depinning 
transition which exhibits prominent hysteresis behaviour 
(Fig. ^ inset) . Such a depinned state is inhomogeneous 

and undergoes plastic flowdl El El El El El 

con- 
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FIG. 4: (a) The power spectrum Ptf>(f) of current fluctua- 
tions in the moving triangle phase (thin line F = 55) and 
in the coexistence region (bold line F = 100), logarithmi- 
cally binned and plotted for a range of frequencies below the 
washboard frequency. Current fluctuations in the coexistence 
region are enhanced, also showing a 1// decay at intermediate 
frequencies, (b) Pn(f) for the fluctuations of the numbers of 
4(D), 5(o) and 6(A) coordinated particles in the coexistence 
region. The number of 7 coordinated particles (not shown) is 
vanishingly small. The straight line in both figures represents 
the 1// law. 



sistcnt with earlier numerical work. For larger F the 
velocity approaches the asymptotic behaviour vcm = F. 

The structure factor S(q) of the plastically moving 
phase (C) obtained just above the depinning transition 
consists of liquid- like isotropic rings (Fig. mi)). Upon in- 
creasing F, the circular ring in S{q) concentrates into six 
smeared peaks which we associate with a hexatic glass 

(D) |lj|. Fig. 121b) shows the evolut ion of the hexatic cor- 
relation function g&{r) — (ip6(0)ipQ(r)) , as F is varied 
across (C) — » (D) — * (E). Note the sharp exponen- 
tial decay of hexatic correlations in (C), the power-law 
or quasi-long-range order (QLRO) decay in (D) and the 
saturation (LRO) of this correlation function in regime 

(E) . The plastic flow regime (C) expands at larger i>3 |l7J • 
On further increasing F, the hexatic glass recrystallizes 
into a structure which depends on the value of v%: for low 
t>3 the final crystal is triangular (E) whereas for large V3 it 
is square (F). For intermediate W3 the system first freezes 
into a triangular structure but subsequently transforms 
into the square via an intervening "coexistence" regime 
(G) best described as a mosaic of dynamically fluctuating 
square and triangular regions. 

The phase diagram of the pure system in thermal equi- 
librium accomodates fluid, triangular solid and square 
solid phases 6]. In the driven system, as F is increased, 
analogous phases appear in approximately the inverse or- 
der to the sequence obtained in the pure case as T is 
increased. This observation agrees roughly with the KV 
proposal [9j , which identifies the shaking temperature T s 
with T. The shaking temperatures are predicted to fall 
as ~ 1/v 2 and as ~ 1/v in the drive and the transverse 
directions respectively, consistent with our observations 



in Fig. Ufa). We find that T" is nearly independent of 
V3. Importantly, within the putative coexistence regime, 
T" behaves non-monotonically, implying a breakdown of 
the KV prediction (see Fig. Ufa)). Typically, for a par- 
ticular disorder configuration and for 5.5 < v$ < 8.5, T S Q 
appears to increase sharply at a well defined F, signify- 
ing the start of coexistence. Within G, T S Q remains high 
but drops sharply at the upper limit of G, to continue 
to follow the interrupted KV behavior. The limits of the 
coexistence region, though sharp for any typical disorder 
realization, vary considerably between realizations. 

Within the coexistence region the probability of ob- 
taining triangular (square) regions appears to decrease 
(increase) roughly linearly with F; see Fig. EJc). Real 
space configurations (Fig. |2J) exhibit islands of square 
and triangular coordination connected by interfacial re- 
gions with predominately 5 coordinated particles. Par- 
ticles with coordination 7 are typically associated with 
dislocations which are scattered randomly in the inter- 
face. This configuration, in the co-moving frame, is ex- 
tremely dynamic, with the islands rapidly interconvert- 
ing between square and triangle. This interconversion 
has complex temporal properties: the power spectrum 
of coordination number fluctuations shows a prominent 
1// fall off over several decades. In addition, particle 
curent fluctuations are enhanced by an order of magni- 
tude, also displaying a regime of 1// behavior (Fig. QJ, 
although over a restricted range as a consequence of the 
proximity to the washboard frequency. This result would 
be hard to interpret in the absence of a genuine coexis- 
tence phase: above the depinning transition, increasing 
the driving force would naively be expected to reduce 
current noise monotonically, as observed in all previous 
simulation work on related models 0, 0, fl4| 

Renormalization group arguments suggest that nei- 
ther translational LRO nor QLRO survive in the dis- 
ordered moving state in two dimensions at any finite 
drive [H El E3, their closest analog being the moving 
Bragg glass state argued to be stable in three dimensions 
and higher j^. The ordered square (F) and triangular 
(E) states we obtain at large drives are then to be un- 
derstood as a finite size effect arising from the restricted 
size of our simulation box, although the crossover length 
scales can be very large at weak disorder^l [2jj. The 
possibility of alternative dynamically stabilized states 
with reduced levels of ordering, such as driven transverse 
smectics, is attractive^! 0]. In contrast to some previ- 
ous we sec no evidence for smcctic order 
and flow in weakly coupled channels at large drives - our 
channels always remain strongly coupled - but note that 
moving states in which channels transverse to the drive 
direction are effectively decoupled may be stabilized at 
higher levels of thermal noise or randomness|14j . 

No general arguments seem to rule out QLRO in the 
"moving hexatic glass" phase in two dimensions and our 
simulations support this possibility; see also Ref. |l6| . 
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The drive-induced stabilization of the triangular lattice 
state we obtain, even well into regimes where V3 would 
favour a square, is an unusual result. Our finding of 
a distinct coexistence regime (G) separate from plastic 
flow states, hexatic and crystal, is novel. Such coex- 
istence occurs in a narrow and reproducible regime in 
parameter space and is characterized by a variety of dy- 
namic anomalies, including enhanced noise signals with 
1// character. The coexistence state appears to be a gen- 
uine non-equilibrium state, separated from other regimes 
through sharp non-equilibrium transitions. 

Theoretical work has, so far, neglected the possibility 
of such dynamic phase coexistence in the non-equilibrium 
steady states of driven disordered crystalsji^- A non- 
disordered but frustrated system closely related to the 
one considered here has been proposed recently as a 
model for the dynamical heterogeniety seen in the glassy 
stately. In this model, fluctuating regions of crystalline 
ordering within a liquid background are argued to be 
responsible for the anomalous dynamical behaviour and 
slow relaxation in the glassy state, a physical picture 
which shares some similarities to our ideas regarding the 
coexistence regime. 

In conclusion, we have shown here that the competition 
between structural phase transitions in a pure system as 
modified by disorder, coupled to the non-equilibrium ef- 
fects of an external drive, can have a variety of non-trivial 
consequences. The ubiquity of structural phase transi- 
tions in the vortex state of a large number of supercon- 
ductors which have been studied recently, as well as the 
relative ease with which the vortex state can be driven, 
suggests experimental situations in which the ideas here 
should find application. Functionalized colloidal particles 
driven over random substrates constitute a novel system 
on which our proposals can be tested 0, . Such sys- 
tems have the further advantage that interparticle inter- 
actions can be tuned both through surface modifications 
and through the application of external fields 23] . 

The authors thank M. Rao, J. Bhattacharya and A. 
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ported by the DST (India). 
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